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SUMMARY 


A theory for the optimum performance of a rotor hovering out of ground 
effect is developed. The performance problem is formulated using general 
momentum theory for an infinitely bladed rotor, and the effect of a finite 
number of blades is estimated. The analysis takes advantage of the fact that 
a simple relation exists between the radial distributions of static pressure 
and angular velocity in the ultimate wake, far downstream of the rotor, since 
the radial velocity vanishes there. This relation permits the establishment 
of an optimum performance criterion in terms of the ultimate wake velocities 
by introducing a small local perturbation of the rotational velocity and 
requiring the resulting ratio of thrust and power changes to be independent 
of the radial location of the perturbation. This analysis fully accounts for 
the changes in static pressure distribution and axial velocity distribution 
throughout the wake as the result of the local perturbation of the rotational 
velocity component. This improvement over earlier theories is shown to have 
finite contributions to the optimum distributions of circulation and inflow 
as well as to the rotor performance. 

The velocities in the ultimate wake are related to the circulation and 
inflow distributions at the rotor disk using a constant contraction ratio for 



the wake. Optimum distributions of the circulation and the inflow, together 
with the performance characteristics are presented for sixteen values of the 
thrust coefficient ranging from 0.001 to 0.050. 


I. INTRODUCTION 


The development of a suitable theory for predicting the optimum performance 
of a rotor hovering out of ground effect (OGE) has long been a central problem 
of rotorcraft aerodynamics. Rotorcraft typically are required to be able to 
hover OGE when occasion demands. Since the power required for a rotor to 
develop a given amount of thrust is the greatest during hover OGE, the payload 
capability of the rotorcraft is often limited by the hovering performance of 
the rotor OGE. For a modern rotorcraft, the payload is likely to be l/5 to 
l/4 of the gross weight at take off. Consequently, a rotor which develops a 
thrust 5$ less than the maximum attainable while hovering OGE will lead to a 
deficiency in payload capability of 20$ to 25$. It is thus of great practical 
importance to study the performance of a rotor hovering OGE. The basic task 
in such a performance study is the determination of a distribution of inflow 
velocity, over a given rotor disk, that leads to a minimum amount of power 
expenditure for a given amount of thrust. This information, when used in 
conjunction with the well known blade-element theory, permits the prediction 
of the performance limit and the required radial distribution of circulation 
in order to approach this limit. It thus provides a rational basis for blade 
design. 

A hovering rotor can be considered as a special case of a propeller, with 
zero advance velocity. Marine and aircraft engineers have long been concerned 
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with propeller aerodynamics . In the majority of situations of interest to 
propeller aerodynamicists , however, the propeller-induced velocity is small 
in comparison with the propeller's velocity of advance. As a result, 
propeller performance theories customarily make use of simplifications that 
are valid only if the induced velocity is small. These theories therefore 
usually do not go over to rotor theories in a straightforward manner and 
rotor performance predictions based on generalized propeller theories do not 
always correlate well with experimental data (Ref. l) . 

In recent years, the search for better methods of predicting rotor 
performance in hover yielded a copious volume of literature on the subject. 
Several recent articles contain comprehensive bibliography (Ref. 2) of the 
earlier work as well as review and reassessment (Ref. 1) of "classical" 
theories. The brief discussion given below provides a summary account of 
previous theories and present efforts. Only a few pertinent articles are 
referred to here. In many cases, a large number of articles exist which are 
based on essentially the same approach. Only one representative article is 
then mentioned. 

The earliest theory on propeller performance was the axial momentum 
theory, based on the actuator disk ( inf ini tely-b laded propeller) concept by 
Rankine (Ref. 3) and Froude (Ref. 4) . Within the context of this theory, the 
propeller operates without any frictional drag on the blade and induces no 
rotation in the slipstream. Since the flow upstream of the propeller disk 
is irrotational, the tangential velocity component at the disk is zero, and 
it follows from the Kutta- Joukowski theorem that the disk cannot be subjected 
to a lift force. Disregarding this inconsistency, two well known conclusions 
were obtained: first, an optimum propeller is one with a uniform axial 
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velocity at the propeller disk, and second, when applied to a hovering rotor, 
the minimum power coefficient is equal to the thrust coefficient to the 3/2 
power divided by the square root of 2. The theory, however, provides no 
information on the radial distribution of circulation over the disk and 
consequently no indication as to how the blades should be designed for optimum 
performance . 

In general, the energy loss due to the rotational motion in the slipstream 
of a propeller is small. The neglect of the slipstream rotation was therefore 
thought to be justified for the purpose of predicting the performance of 
propellers. The actual performance of propellers, however, generally does 
not meet the expectation of the axial momentum theory even when the various 
loss factors, not including the slipstream rotation, are taken into account. 
This deficiency is explained by the general momentum theory, which incorpo- 
rates a procedure to account for some of the effects of slipstream rotation. 
One of the major conclusions of the general momentum theory is that the 
slipstream rotation, though small in terms of the energy content, exerts an 
important influence on the optimum distribution of inflow velocity over the 
propeller disk. A detailed presentation of the general momentum theory is 
given in Ref. 5 which also contains a presentation of efforts directed toward 
the study of helicopter rotors. The general momentum theory does provide 
information on the optimum distribution of circulation over the disk. The 
solutions as given in Ref. 5, however, are approximate and only partially 
account for the effects of slipstream rotation. 

The theory which represents the current state of the art for routine 
calculation of propeller performance is the vortex theory. The basic vortex 
theory represents the wake trailing the blades by a distribution' of concentric 
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cylindrical vortex sheets. These vortex sheets describe the radial variation 
of the blade circulation and thus imply an infinite number of blades . The 
effect of a finite number of blades is obtained approximately by Prandtl 
(Ref. 6) using a tip loss factor. An improved analysis by Goldstein (Ref. 7) 
represents the trailing vortex sheets, for optimum performance, by discrete 
helicoidal surfaces of constant helix angle and moving as rigid surfaces. 
Goldstein's work formed the basis of much of the subsequent analyses by Lock, 
Theodorsen, Lerbs, etc., (Refs. 8, 9> 10, and 11) that proved to be suffi- 
ciently accurate for predicting propeller performance. The use of the vortex 
theory for predicting hovering rotor performance, however, yielded overly 
optimistic results (Ref. 1) . Reviews of the vortex theory suggest that the 
wake contraction, the slipstream rotation, and the associated non-uniform 
inflow are the factors that contribute most to the inability of the vortex 
theory to predict the hovering rotor performance accurately. This contention 
is supported by recent results obtained from the numerical computation of the 
entire rotor-induced flow field. Before discussing the numerical approach, 
it is to be noted that the work of Lerbs (Ref. 11) extends the vortex theory 
to heavily- loaded propellers where the effect of wake contraction is important. 
Lerbs' work, though well known in marine engineering, has not been extended 
and applied to the study of rotor performance. As is the case of the general- 
momentum theory, the work of Lerbs is based on certain simplifying assumptions 
that cannot be justified for a hovering rotor. The modification of Lerbs' 
analysis for rotor applications is by no means straightforward. 

The continuing requirements of higher forward flight speed and larger 
rotorcraft, coupled with the need to keep the rotor reasonably small from 
weight and operational flexibility considerations, led to the development of 
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rotors with higher rotational speed and disk loading. This resulted in 
increased inaccuracy of the existing theories in predicting the hovering 
performance. The most recent efforts to overcome this difficulty is the 
development of numerical methods for the integration of the Biot-Savart law; 
the use of which permits the computation of the entire flow field induced hy 
the rotor. It is noted that calculation of the rotor performance requires 
only a knowledge of the inflow velocity over the propeller disk. The quadra- 
ture of the Biot-Savart law, in order to establish the inflow velocity 
components at the disk, however, requires the computation of the entire vortex 
system at and trailing the disk. Since the computation requires an iterative 
procedure and the number of data points involved in each iteration is large, 
the approach is time-consuming even by modern computational standards . To 
solve a problem where the propeller geometry, rather than the blade circula- 
tion distribution, is specified, an additional iterative procedure is needed 
to establish the circulation distribution corresponding to the specified 
propeller geometry. Further, to obtain the propeller geometry that would 
give the optimum performance, a parametric study involving a large number of 
geometries must be treated. For this reason, although some success has been 
reported in the development of the Biot-Savart law approach, the utility of 
this approach for design purposes is limited at the present. 

The purpose of this report is twofold, (a) to present a simple theory, 
together with numerical results , for rotor performance and (b) to clarify 
certain features of the performance problem pertinent to rotors hovering OGE. 
The theory presented here is similar to the general momentum theory described 
in Ref. 5 in that the rotor is represented by an actuator disk. It is there- 
fore subjected to some of the limitations of the general momentum theory. In 
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particular, the correction factor for finite number of blades must be 
estimated separately, for example, by an extension of Lerbs’ (Ref. 10) method 
for heavily loaded propellers. The present method, however, is more complete 
in that it fully accounts for the effects of slipstream rotation. This 
improvement, which is not important in the case of a lightly loaded propeller, 
has finite contributions to thrust and power requirements of a rotor. The 
use of the present theory permits the calculation, in a straightforward 
manner, of the distribution of circulation and inflow velocity over the rotor 
disk for the infinitely bladed case. These distributions are, along with the 
figure of merit and the optimum power coefficient, presented for several 
values of the thrust coefficient. 
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SYMBOLS 


a = u^/OR^, dimensionless axial velocity in the ultimate wake 

A value of "a" at the rim of the ultimate wake 

b = dimensionless angular velocity in the ultimate wake 

B value of 'V at the rim of the ultimate wake 

H total head of the fluid 

H total head of the ambient fluid 

o 

H-^ total head of the fluid in the ultimate wake 

k = r^/r, local contraction ratio 
K = R-j/R, overall contraction ratio 
M figure of merit, Eq. (51) 

n a constant in the optimum performance criterion, Eq, ( 32 ) 

R = n/R^, a dimensionless constant 
R^ number of blades 

■** 

p static pressure of the fluid immediately upstream of the rotor disk 

P o static pressure of the ambient fluid 

p^ static pressure of the fluid in the ultimate wake 
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p 7 static pressure jump across the rotor disk 

P power expended by the rotor 

P^ power coefficient, Eq. (50 ) 

Q torque on the rotor 

r radial position on the rotor disk 

r^ radial position in the -ultimate wake 

R radius of the rotor 

R^ radius of the ultimate wake 

S area of the rotor disk 

area of the ultimate wake 

T thrust developed by the rotor 

T thrust coefficient, Eq. (49) 

c 

u axial velocity at the disk 

axial velocity in the ultimate wake 
x = r-j/R-^, dimensionless radial position in the ultimate wake 
y = r/R, dimensionless radial position on the disk 
T circulation at the rotor disk 
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circulation of each blade 


an gular velocity of the fluid immediately downstream of the rotor 
disk 

angular velocity of the fluid in the ultimate wake 
angular velocity of the rotor 


density of the fluid 



II. FORMULATION OF THE HOVERING ROTOR PROBLEM 

The equations governing the performance of hovering rotors are derived 
below within the context of the general momentum theory. The general flow 
features are shown in Fig. la,. The rotor is represented by an actuator disk 
of radius R and driven to rotate about its axis at a constant angular velocity 
Q in an incompressible, inviscid fluid. The flow is steady and symmetric 
about the axis of the actuator disk. The rotational component of the velocity 
is, however, non-zero in the slipstream. The fluid outside the slipstream is 
irrotational and is at rest at infinity. The total head of the fluid outside 
the slipstream is therefore equal to the static pressure of the fluid far 
upstream and is a constant. The actuator disk imparts an increase in total 
head to the fluid passing through the disk in the form of an abrupt static 
pressure increase and a change in the angular velocity of the fluid from the 
zero angular momentum value immediately upstream to some finite value 
immediately downstream of the disk. The law of mass conservation requires 
the axial velocity component to be continuous across the disk. The radial 
velocity component is taken to be continuous across the disk within the 
context of the general momentum theory. 

The slipstream, in which the flow is rotational and possesses a higher 
total head than the flow outside, contracts downstream of the disk and forms 
an ultimate wake far downstream of the disk. In the ultimate wake, the 
radial velocity component vanishes and the flow properties are independent 
of the axial coordinate. Let the flow through the annular element dS = 2nrdr 
at radius r at the disk, shown in Fig. la pass through the annular element 
dS-^ = Pnr^dr^ in the ultimate wake. The subscripts "1" designate flow 
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conditions in the ultimate wake. The unsub scripted variables designate the 
corresponding flow conditions immediately downstream of the disk. 

The law of mass conservation states that 

u^r^cfr^ = urdr (1) 


The tangential component of the equation of motion requires the angular 
momentum of the fluid to be constant along streamtubes in the slipstream. 
Thus one has 


(u r 
11 


= tor 


( 2 ) 


The radial- component of the equation of motion requires the radial- 
pressure gradient in the ultimate wake to be balanced by the centrifugal- 
force on the fluid, since the radial velocity component is zero in the 
ultimate wake . Thus 


dp 2 

dr-^ 1 1 


(3) 


Bernoulli’s equation, applied to the flow upstream of the disk gives 


H 


o 



= P 


,i,2 2 N 

+ 2 P (U + V ) 


( 4 ) 


% 

where the subscripts "o" designate flow conditions far upstream of the disk 
and p is the static pressure of fluid immediately upstream of the disk. 

The total head immediate downstream of the disk is 


H = p 


n n 

+ p + -g-p (u + v + 


u) 2 r 2 ) 


(5) 
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where p ' is the pressure jump across the disk and tu is the angular velocity 
of the fluid immediately downstream of the disk. Applying Bernoulli’s 
equation to the flow downstream of the disk gives 

= p + p 7 + ip (u 2 + v 2 + cu 2 r 2 ) 

= P x + ip (u/ + u^ 2 ^ 2 ) (6) 

Equation (5) and (6) gives 

p' - (% - H q ) - ipm 2 r 2 (7) 

The element of torque of the disk is equal to the flux of angular momentum 
imparted to the fluid passing through the annular disk element dS. Thus 

dQ, = 2rrpuu)r^dr (8) 

The element of power expended hy the disk therefore is 

dP = HdQ 

= 2rrpQuu)r^dr (9) 

This element of power is also equal to the increase in total head per unit 

time imparted to the fluid passing through the annular disk element, i.e., 

s 

dP = SIttCe^ - H ) urdr 


Consequently 



pQur = pQtu^r^ c 


( 10 ) 
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Putting Eq. (10) into Eq. (7) yields 

p ' = p(0 - -gts) u)r 2 (11) 

The element of thrust developed by the, disk is 

dT = p 'dS 

= 2TTp(n - |<jo) u)r^dr (12) 

The total thrust developed and the total power expended hy the entire disk 
are expressible in terms of the velocities at the disk by using Eqs. (9) and 
( 12 ) : 

( 13 ) 

o 

rR o 

and T = 2rrp (n - g-cu) uor^dr (14) 

o 

In terms of the velocities in the ultimate wake, the total power is, 
using Eqs. (1), (2), and (13 ), 

P = 2npn J^ 1 u i (jU 1 r 1 ^ dr 1 (15) 

o 

* 

An expression for the total thrust in terms of the wake properties may be 
obtained by specializing the expression for a propeller in forward motion 
given in Ref. 5 to the case of zero advance velocity: 


= 2npQ 


3 

uu)r dr 
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( 16 ) 


I = 2n[p f 1 u 1 2 i- 1 dr 1 - J® 1 (p - Pj) r x dr 
o o 

Although the derivation presented in Ref. 5 requires certain assumptions not 
appropriate for the hovering rotor case, it is shown in Appendix A of this 
report that Eq. (16) is valid for the hovering rotor. 

Using Eqs. (4) , (6) , and (10) , one obtains 

P Q - Pq = ipiiq 2 - p(fi - iwqjajqr ^ 2 (IT) 

Equation (16) can therefore be re-written as 

T = 2n p 1 [?u 1 2 + (o - fto-jj u> 1 r 1 2 ] r 1 dr 1 (18) 

o 

Differentiating Eq. (17) with respect to r^ gives, after using Eq. ( 3 ) 
and re-arranging terms 



At the rim of the ultimate wa%e, r^ = R^ and p Q = p^. Equation ( 17 ) thus 
gives 

u 1 2 (Ri) = [20 - ^(Rq)] ^(R^Rq 2 (20) 
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III. CRITERION FOR OPTIMUM PERFORMANCE 


In the ultimate wake where the radial velocity vanishes, the simple 
relationship between the pressure and the angular velocity, Eq. (3) , makes 
it possible to express the axial velocity Up as a function of the angular 
velocity 0)^; i.e., if U)p(rp) is 
quadrature of Eq. (19) • Such a simple relationship does not exist at the 
disk, where the radial velocity component is non-zero. Consequently, in 
deriving a criterion for optimum performance, it has not been possible to 
consider the flow conditions at the disk directly. In the following analysis, 
a criterion for optimum performance is established in terms of the velocity 
distributions in the wake. These velocity distributions will subsequently 
be related to the velocity distributions at the rotor disk and to the thrust 
and power coefficients of the rotor. 

Consider a given distribution of angular velocity cu^(r^) . Let this distri- 
bution be perturbed to U) p( r p) + Aujp, with Au)p given by 

'0 , 0 < r x < C 

AcUp = l e , £ < t ± < G + AC (21) 

1 0 , c + AC < rp < Rp 

where e « cOp(G) and AG « £. As & a result of this perturbation in u)p, there 
exist perturbations in the pressure and the axial velocity in the ultimate 
wake as shown in Fig. lb. From Eq. (3) , one has 

1 p«Bp 2 5cL? (22) 

1 

Thus, to the first order in AWp/oUp, the perturbation in pp, which results 



specified, then u p( r p) follows by the 
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from the perturbation in ou^, is 


AP 1 = - 2p J^o^doijgd? 


(23) 


Placing Eq. (21) in Eq. (23) gives 


f - 2po) 1 (OC(eAO 


A Pl = 


0 < r ± < C 


- 2 P (ix L (0C[e(C + AC - r x )] C < ^ < C + AC 


C + AC < r 1 < R x 


(24) 


where only the lowest order term in e/au^ and in AC/C are kept in each interval 
of interest. 

Using Eq. (17) , the perturbation in u^ can be expressed as 


2 

pu-^Au-^ = - Ap 1 + p(n - r 1 Ac « 1 


(25) 


Thus from Eqs. (21) and (24) 


' 2p c» 1 ( C) C(®AC) 

pu-j^A^ = ^ 2pci) 1 (C) C[ e (C+AC-r 1 ) ] + p[n-(D 1 (C) ]c 2 e 
0 


0 < r 1 < C 


C < r 1 < C + AC (26) 


C + AC < r ± < R 1 


Using Eqs. (16) and (15) , the changes in thrust and in power due to the 

% 

perturbation in o>^ are 


<»R 


AT = 2 rr 


1 " 


2p u^Au^ + 


(27) 
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and 


(*1 o 

AP = 2npfi J (^Au^ + u> Au ) r^di^ (28) 

o 

Using Eqs. (21), (24), and (26), one obtains 

AT = 2rrp[2n - <^(0 ] £ 3 (eA£) (29) 

and 

d§} (eA£ 

Equations (29) and (30) may be combined to express AP in the form 

AP = g(Q AT (31) 

where g(£) is a function of £ , the location where the angular velocity 
perturbation is introduced. Suppose that an angular velocity perturbation 
is introduced at £^, with resulting changes in power and thrust AP^ and AT-^. 
Suppose a second angular velocity perturbation is introduced at £g with 
resulting changes in power and thrust APg and ATg. By suitably selecting 
the values of eA £ at £^ and £g, the net change in thrust due to the two 
perturbations is made to zero, i.®., AT = AT^ + ATg = 0. Then, using Eq. ( 31 ) , 
the net change in power is 

AP = AP 1 + APg = [gty - g ( £g) ] AT x 

Since AT-^ can be made either positive or negative by choosing the sign of ”e" 



r o ^(£)£ 5 (n-%) 

AP = 2TT P n |u 1 £- 3 + + 2U) 3 


'o U 1 


18 



in Eq. (29) , it is possible to obtain a negative value for AP if g(C^) f g(Cg) • 
In other words, it will be possible for the rotor (with fixed R and 0) to 
develop a given amount of thrust while expending a smaller amount of power by 
altering the given angular velocity distribution u>^(r.j) . Ehis distribution 
therefore is not optimum. 

It is obvious, therefore, that the optimum distribution ux^(r^) requires 
g(C 1 ) = g(£g) • •^ n °ther words, the criterion for optimum performance is that 

g(£) is a constant. Let this constant be On. From Eqs. (29) and (30) , the 
optimum performance criterion is 


? 2 f^l q II 

“i (2n - “1) * ¥1 + 2u) i J u 7 5 + ST r i ^ 


(32) 


A differential form of this performance criterion is obtained by dif- 

ferenting both sides of Eq. (32) with respect to r^. !Ehe resulting equation 

and Eq. (19) > both containing the derivatives of u^ and o>^ with respect to r^, 

du^ d(JD^ 

permit each of the derivatives, and rgy- , to be expressed as a function of 
u^, u)^, and r^. In terms of dimensionless variables x = r-J R^, a = u^/QR^, 
b = u^/Q, and the dimensionless constant N = n/R^, the two derivatives are 


da 

dx 


- 2ab(l - b) [2a 2 - a(4 - b) W + 2bx 2 ] x 
a^ - 2a 3 N - a 2 b(l - 2b) x 2 + b 2 (l-b) 2 x 4 


(33) 


and 


db _ 2b[a^ - a 3 (2-b) N + a 2 b(3 - 2b) x 2 - b 2 (l - b) 

^ [a 4 - 2a% - a 2 b(l - 2b) x 2 + b 2 (l - b) 2 x 4 ] x 


(34) 


At the rim of the slipstream, a and b are related by Eq. (20) which, in the 


19 



non-dimensional form, is 


A 2 = B(2 - B) (35) 

where A and B are values of a and b respectively at x = 1. 

Taking the upper limit of integration to be R^ in Eq. (32) and using 
Eq. ( 35)5 one obtains an expression for N: 

< 36 ) 

O 

For any given value of B (or A) , the system of equations (33) to (36) 
determine functions a(x) and b(x) that give optimum rotor performance as 
well as the value of N. For a given value of B, the value of A is determined 
by Eq. (35) • The value of N depends on the functions a(x) and b(x) and is 
determined in an iterative procedure together with the functions a(x) and 
b(x) . For the first iteration, an estimated value of E is placed in Eqs. (33) 
and (3^) , which are solved for a(x) and b(x) using a fourth order Runge-Kutta 
method for simultaneous first order ordinary differential equations . The 
calculation begins at the rim, x =1, where the boundary values A and B are 
known, and proceeds inward towards the axis, x = 0. The results are used to 
obtain a corrected value of N for the subsequent iteration. The iterative 
procedure continues until the quadrature of Eq. (36) yields a value of N 
differing insignificantly from the input value of N in Eqs. (33) and (3V) . 

The detailed computational procedures are presented in Appendix B. 
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IV. PERFORMANCE AND FLOW CONDITIONS AT THE ROTOR DISK 


For a given radius of the ultimate wake, R^, and a given rotor angular 
velocity 0 , the optimum performance functions a(x) and b(x) specify the 
optimum flow field in the ultimate wake. The corresponding power expended 
and thrust developed by the rotor are determinant from Eqs. (15) and (18) . 

The optimum circulation and the axial inflow velocity at the disk as well as 
the power and thrust coefficients and the figure of merit of the rotor can be 
computed based on the known functions a(x) and b(x) , provided that a func- 
tional relationship between r and r^ can be established. In theories of 
rotor and propeller performance, it is customary to obtain a differential 
form of the thrust element by differentiating equation ( 18 ) : 

“ = 2tt p [* u/ + (0 - i«. x ) Vl 2 ] (37) 

Equating the right sides of Eqs. (12) and (37) gives 

dr [ u l 2 + ( 2 " - "l) Vl 3 ] r l 

= o (38) 

ar l (20 - m) U)r^ 

Equations (2) and ( 38 ) , together with the boundary condition 

r = 0 at r, = 0 ( 39 ) 

*. 1 

determine a functional relationship between r and r^, including the overall 
contraction ratio K = R-^/R, when the velocity distributions u^ and in the 
wake are prescribed. 

Equation (37) obviously is not the only expression that satisfies Eq. (18) ; 
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for any expression g(r 1 )dr 

o 

right side of Eq. (37) , and the resulting expression still satisfies Eq. (18) . 
In fact, it is shown in Appendix C that for bounded and r(r-j) a one to one 
function, no solution of Eq. (38) exists which satisfies the boundary condi- 
tion (39) • A more detailed discussion is given in Ref. 15. 

The function r(r^) can be obtained as a part of the solution of the 
complete flow field between the disk and the ultimate wake. The differential 
equation governing the flow field is elliptic (Ref. 12) and it has been 
suggested (Ref. 13) that the exact contraction ratio can be obtained only by 
the complete solution of the governing equation. Such a complete solution 
constitutes a major task demanding a large amount of computational effort. 
Accordingly, in the present work, an approach which requires neither the use 
of Eq. (38) nor the solution for the entire flow field is developed. In this 
approach, the local contraction ratio k = T-Jr is taken to be independent of 
the radial position r (or r-j) . 

It is recognized that in general the local contraction ratio depends 
on the radial position. The results of Ref. 10, however, indicate that for 
heavily loaded free-running propellers, the contraction ratio is nearly 
independent of the radial position, except near the axis of the propeller, 
where the dependence of the contraction ratio on the radial position has only 
small effects on the propeller ’s performance. Thus it is expected that the 
assumption of a constant contraction ratio is reasonable for rotor performance 
studies . 

With k(r) = K, one has 

r = r^/K and dr = dr^/K ' (40) 


^ which gives J g(r-^) cLru = 0 can be added to the 
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Thus Eqs. (l) and (2) become 


and 


v 2 

u = u^K 
v 2 

w = od-j^K 


(41) 

(42) 


The value of K can be determined by equating the right sides of the two total 
thrust equations, (14) and (18) . Using Eqs. (40) to (42) and rearranging 
terms, one obtains the following expression for K: 


K = 


fs : 


r l ar l 


1 + 


20 f 1 o) 1 r 1 3 dr 1 

o 


(43) 


Integrating by parts gives 


rR 


1 2 , U ,‘4s,> K / J * 

U 1 r l dr l “ 2 2 


1 2 du l - 

r l d?“ ^1 
o 1 


(44) 


Using Eqs. (19) and (20) and integrating by parts yields 


Jt 


^ u^r^dr^ = 2n 


r*R- 


1 3. 

u^r^ dr^ 


o 


o 


Putting Eq. (45) into Eq. (43) gives 

K = l//2~ 


(45) 


(46) 


which is the value predicted by the axial momentum theorem. 

The axial inflow velocity and the circulation at the disk are related to 
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the axial and angular velocities in the ultimate wake by 


and 


u(r) = i u-^r-j) 

r(r) = 2rrw^ ( r^) 


The thrust and power coefficients are given by 


T 

c 


T 

ttr 2 p (Or) 2 



2 

a 


+ (2 - b) 



P 

c 


P 

ttr 2 p (Or) 3 


1 

2J2 


rl 


abx 3 dx 


The figure of merit is given by 


T 3//2 

M = — 2 

,/F P 

C 


V. RESULTS MD DISCUSSIONS 


( 47 ) 

( 48 ) 


( 49 ) 

(50) 


( 51 ) 


The optimum performance characteristics of the hovering rotor are 
presented in Table 1. The minimum power coefficient, the optimum figure of 
merit and the non-dimensional circulation at the disk edge are given for 16 
values of the thrust coefficient between 0.001 and 0.050. This range of 
thrust coefficient was considered to bracket the normal operating range of 
hovering rotors. 

The optimum radial distributions of circulation on the disk, are presented 
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in Fig. 2 for the 16 values of the thrust coefficient. It can he seen that 
for small values of T , the optimum circulation is nearly independent of the 
radial position except near the axis. For higher values of T , the deviation 
of the optimum circulation from a constant value extends over a larger region 
and the magnitude of the deviation is greater. 

The corresponding optimum distributions of the axial inflow velocity at 
the disk are presented in Fig. 3* It can be seen that a similar pattern 
exists for the axial inflow velocity as for the circulation. That is, the 
optimum axial inflow velocity is nearly uniform for small values of T c except 
near the axis. The deviation of optimum axial inflow velocity from a 
constant value is more significant for the higher T values. The slope of 
the curves u(r) is zero at r = 0 as expected. 

In Ref. 5s an approximate solution for the hovering rotor problem is 
given based on the assumptions that the radial pressure gradient in the 
ultimate wake is negligible and that the non-dimensional velocities at the 
disk u(y) = ^ and u>(y) = q are 




'V 

u 


_5SL 


2 2 
m + y 


(52) 


and = — 2 ~— — (53) 

m + y 

where m is a constant and y = r/R. The approximate function cu(y) and u(y) 
are compared to the present results for optimum performance in Figs. 4 and 5 
for the case T c = 0.01. A similar comparison for the circulation is given in 
Fig. 6. The comparisons show that the approximate results deviate signifi- 
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cantly from the present results, particularly in the inboard region of the 
disk. The deviations are more severe for larger values of T . The approxi- 
mate results give excessively high value of cu and zero value of u at the axis, 
as is shown to he incorrect in Appendix C. The present solution gives a 
moderate value of w and a finite value of u at the axis. 

The variation of pressure in the ultimate wake is plotted against the 
radial position in Fig. 7> again for the case T c = 0.01. As can be seen from 
Eq. (17) , if the pressure gradient in the ultimate wake is neglected, as was 
done in Ref. 5? the axial velocity vanishes at the axis. If the pressure 
gradient is taken into account, as is in the present analysis, the axial 

p 

velocity is finite at x = 0. The value of (p - p^) /p (fiR^) = 0.0123 at 
x 1 = 0, given in Fig. 7, corresponds to a value of u^ = O.I 57 Qr^, which is 
comparable to the maximum axial velocity in the wake, 0.207 fiRq. The effect 
of the pressure gradient in the ultimate wake on the thrust developed by the 
rotor is represented by the second term on the right side of Eq. (16) . The 
neglect of the pressure gradient would lead to, according to Fig. 7> a 
prediction of the thrust coefficient approximately 2 % higher than the optimum 
for the case T c = 0.01. More important is the fact that the optimum radial 
distributions of circulation and axial inflow velocity predicted by the 
approximate theory, which neglects the pressure gradient, contain significant 
errors . 

In Fig. 8, the optimum power coefficient is presented as a function of the 
thrust coefficient. The dependence of the figure of merit on the thrust co- 
efficient is shown in Fig. 9* The present result is compared with the approxi- 
mate result of Ref. 5 and the ideal figure of merit. The figure shows that 
the approximate method of Ref. 5 predicts the rotor performance somewhat more 


26 



optimistically than the present theory. Both the approximate method and the 
present method predict performance substantially lower than ideal. The 
deviation from the ideal performance increases with increasing thrust 
coefficient . 

The present results are primarily concerned with an infinitely bladed 

rotor. The effect of a finite number of blades may be determined, for example, 

by extending herbs' work (Ref. 10) on heavily loaded, free running propellers 

to hovering rotors. An estimate of the effect of number of blades on the 

optimum distribution of circulation is provided by modifying Prandtl's 

approximate method (Ref. 6) of calculating the "tip loss" factor and applying 

the result to the hovering rotor case. The procedure is outlined in Appendix 

D. The results are shown in Figure 10 for rotors with 2, 4, and 6 blades at 

T = 0.010. Based on these results, the corrected thrust coefficients for 
c 

2, 4, and 6 blades are respectively 0.00759? 0.00865? and .00906. The 
corresponding corrected power coefficients are .000547? 0.000624, and 
0. 000655 • The figures of merit are 0.855? 0.9H? and 0.931* In comparison, 
without the tip loss correction, the power coefficients for the infinitely 
bladed rotor are respectively (at T c = 0.00759? O.OO 865 ? and 0.00906) 

0.000479? O.OOO 585 ? and 0.000627. The figures of merits are 0.975? 0.973? 
and 0.972. 
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VI. CONCLUSIONS 


A theory for the optimum performance of a rotor hovering out of ground 
effect is presented. The formulation of the problem is based on the general 
momentum theory of an actuator disk (infinitely bladed rotor) and makes use 
of the fact that in the ultimate wake, far downstream of the rotor disk, the 
radial component of the velocity vanishes and there exists a simple relation 
between the radial distribution of static pressure and that of rotational 
velocity component. The present analysis, which leads to an optimum perform- 
ance criterion, is more complete than previous analyses in that the present 
work fully accounts for (l) the effect of slipstream rotation and (2) the 
existence of a radial pressure gradient in the ultimate wake. 

Numerical results are presented for the radial distributions of circula- 
tion and the axial inflow velocity at the disk, the optimum power coefficient 
and the optimum figure of merit for sixteen values of the thrust coefficient 
covering the normal range of interest in rotorcraft application. Comparisons 
are made between previous approximate results and the present more exact 
results . 

“Several conclusions of the present study are summarized below: 

1. The customarily accepted expression for the differential element of 
thrust in terms of the flow conditions in the ultimate wake, used in previous 
analyses to relate the wake flow to the flow at the disk, is inconsistent with 
the fact that a radial pressure gradient exists in the ultimate wake. 

2. The existence of the radial pressure gradient in the ultimate wake 
has a significant effect on the optimum radial distribution of circulation 
and the optimum performance of the rotor in the range of thrust coefficient of 


28 



interest to rotor craft application. 

3. Numerical results obtained based on the present theory indicate that 
previous approximate theories yielded overly optimistic results on the per- 
formance of hovering rotors. The inaccuracy in the approximate results 
increases with increasing thrust coefficient. As has been noted by earlier 
investigators comparing experimental data with theoretical results, the 
existing theories generally overestimate the hovering rotor performance in 
the higher thrust coefficient range. Thus the correction factor based on the 
present theory may be utilized to improve the previous predictions. 

4 . With the assumption of uniform contraction, the contraction ratio is 
l//2~ , the value predicted by the axial momentum theory. 

5 . The approach developed is potentially useful for predicting the 
optimum performance and the related distribution of circulation of heavily 
loaded propellers and rotors in very general types of axi-symmetric flow, 
such as the shrouded propeller and the ducted fan flows . The effect of finite 
number of blades can be determined separately by modifying available theories. 

•# 
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TABLE 1 


OPTIMUM PERFORMANCE CHARACTERISTICS 


2 

1 x 10 

p x 10 3 

M 

B x 10' 

c 

c 



0.1 

.02246 

.9954 

.2026 

0.2 

•06377 

.9917 

.4091 

0-3 

.1175 

.9885 

.6189 

0.4 

.1815 

.9854 

.8317 

0.5 

.2545 

.9824 

1-0473 

o .6 

• 3355 

.9796 

1.2657 

0.7 

.4239 

.9769 

1.4866 

0.8 

.5193 

.9743 

1.710 

0.9 

.6213 

.9717 

1.9358 

1.0 

.7296 

.9692 

2.1640 

1.5 

1.3569 

.9574 

3-3393 

2.0 

2.1134 

• 9463 

4.569 

2-5 

2.9865 

• 9359 

5.8510 

3-0 

3-9682 

.9259 

7.1841 

4.0 

6.2379 

.9086 

10 .0024 

5.0 

8.8958 

.8887 

13.0246 
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APPENDIX A 


TOTAL THRUST IN TERMS OF THE WAKE PROPERTIES 

Consider a spherical control volume of radius ft centered at the center of 
the rotor disk. Let this control surface he divided into two parts, S and 

SI 

S b , S & being the part of the control surface where the slipstream leaves the 
control volume and S^ being the remaining part of the control surface. The 
momentum theorem gives, for the present steady flow problem, 


' 2— — r 2 _ - 

= ^ (p + pu ) n-i dS & + j (p + pu ) n*i dS^ 


(A-l) 




where p and u are the static pressure and the axial velocity of the fluid, 
n is the unit outward normal and i is the unit vector in the axial direction. 
At large distances from the rotor disk, the effect of the rotor on the fluid 
outside the slipstream is equivalent to that of a sink located at the center 
of the rotor disk. The strength of the sink is equal to the volume rate of 
flow passing through the disk and is finite. In the limit as ft -* ro , on 

# Q 

the magnitude w of the velocity vector goes to zero; i.e., w -» — x where C is 

4 5 % 

a constant. Consequently (p - p) = § pw^ -*» on S, as ft -* ®. The second 

° 2ft 4 b 

integral in Eq. (A-l) thus gives -p (tir^) as ft -* 00 , and equation (A-l) 
becomes 


& 


T = 2tt 


1 [(Pp - P ) + PUp 2 ] r 1 dr 1 


(A- 2) 
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APPENDIX B 


COMPUTATIONAL PROCEDURES 


Equations (33) and (34) constitute a system of two first order ordinary 
differential equations for the functions a(x, N) and b (x, N) . Equation ( 36 ) 
is an auxiliary equation relating N to the functions a(x, N) and b(x, N) . 
Equation (35) relates the values of a and b at x = 1. Thus, if either 
A = a(l) or B = b(l) is given, the functions a and b are obtainable from 
Eqs. (33) 5 (34) , and ( 36 ) by using an iterative method to establish the value 
of N. 


For a selected value of B, the value of A is calculated using Eq. (35) • 

The procedure for the ith iteration is as follows : 

1. With a selected value of , designated , Eqs. ( 33 ) and (34) 

are solved for (x) and b^ (x) using a fourth order Runge-Kutta method 

for simultaneous solution of the differential equations. The computation 

begins at x = 1 where the values of a and b are known and proceeds inboard, 

at intervals of Ax = 0.01, to x = 0. 

* 2 . With the calculated a^ (x) and b^ (x) , a new value of N ^ , 

designated N , is computed by numerical quadrature of the integral in 

Eq, ( 36 ) using a Newton-Cotes formula. 

3 . The difference F^ = N ^ - N. ^ is calculated and compared to 

out m 

EL . If the magnitude of F^ is less than 10 ^ N^^ 5 then the computation 

(l) (l) 

is terminated and the functions a^ ' and b v ‘ are accepted as the optimum 


function a(x) and b(x) for the given value of B. Otherwise a value of N 
is determined for the subsequent iteration. 


(i+1) 


m 
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for the initial iteration is critical. 


The selection of the value N. ^ 

in 

As is shown in Fig. 3? the function a(x) is positive for optimum performance. 
If the value is improperly selected, the calculated function a^ (x) 

may vanish at some points in the interval 0 < x < 1. The numerical quadrature 
of the integral in Eq. ( 36 ) then becomes inaccurate. In actual computation, 
it was found that a value of only slightly smaller than the correct 

value of N may lead to computational difficulties . 

A reasonable value of is obtained by using the form of a and b 

given by Glauert’s approximate expressions, Eqs. ( 52 ) and (53)* That is, 

C,mx 2 

a = 2 (B' 1 ) 

m + x 

„ 2 
Cm 

and b = ^ (B-2) 

m + x 

To satisfy Eq. (35) for all values of m, one has C-^ = = 2. Placing (B-l) 

and (B-2) in Eq. ( 36 ) gives 


N. 


(1) V B ( 2-B) (3-B) 


m 


(2-B ) 1 


(B-3) 


It was found that a linear iteration procedure, i.e., letting 
i +1 i 

leads to very slow convergence. In the actual computation, for 
the second iteration, the value of If. ^ is taken to be O .98 H. ^ . For all 
subsequent iterations, N. is obtained using the secant method, i.e., 
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b(0) = 2 - | (B-6) 

Equations (B-5) and (B-6) provide convenient means for checking the 
accuracy of the solutions. 

In order to obtain the functions a(x) and b(x) for given values of the 
thrust coefficient, instead of given values of B, it is necessary to iterate 
for different values of B. This iteration is straightforward and involves 
no computational difficulty. 
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APPENDIX C 


RESTRICTIONS DUE TO RADIAL PRESSURE GRADIENT IN TEE WAKE 


Equations (19) and (20) give the following equation: 


2 2 2 

u.^ = (20 - cD-j^) + 2 J r l dr l 


(C-l) 


Since iu^ is bounded and non- negative, the first term on the right side of 
Eq. (C-l) approaches zero as r^ -* 0. Thus u-^ is non-zero at the axis = 0. 
The integral in Eq. (C-l) is positive and increases as r, -*0. Hence there 
is a point S and a positive constant K such that for r^ < S, 


u-j 2 + (20 - uo 1 ) u^r 2 2= K 


(C-2) 


Consequently, from Eq. ( 38 ) , one has for r, < S 


dr 

dr 


Kr, 


1 (20 - 10 ) aur^ 


(C~3) 


Sinpe r(r^) is a one to one function, neither uu nor (20 - u)) can be negative 
for r, ^ S. Equation (C-3) may therefore be rewritten as 


, Kr, 
dr ^ 1 


Kt 


K 


dr l ate- 3 " afcy/r 


(o-4) 


where M is a constant such that u>^ < M for r^ < S. Let r be given by 
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K 



2fiMr-jjr ‘ 



(C-5) 


and r / = r-^S) at = S (C-6) 

then, from Eq. (C-4) , one has 

r f ^ r for < S (C-7) 


However, Eqs. (C-5) and (C-6) give 

r /2 = ri 2 (S) - ^ ^(s/rj (C-8) 

r QMr 2 (S)-, 

Consequently r' =0 at ^ = S exp|^ j? J > 0. Since r ^ r' for r 1 < S, 

r ^Mr 1 2 (S) 1 

one has r = 0 at a point r^ ^ S expl ^ > 0 . The condition r = 0 

at r^ = 0 therefore cannot be satisfied. 

It is noted that, if the radial pressure gradient in the ultimate wake is 
neglected, as was done in Ref. 5 j then the last term in Eq. (C-l) vanishes. 
Consequently, u^ = 0 at r^ = 0 and the condition r = 0 at r^ = 0 can be 
satisfied. The consideration of the effects of the slipstream rotation, 
however, requires the pressure gradient term to be retained. 
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APPENDIX D 


ESTIMATED TIP LOSSES BASED ON PRANDTL'S METHOD 

Prandtl's approximate formula for the tip loss correction factor for the 
blade circulation is 

F * | cos" 1 (e" f ) (D-l) 

Tfcr 

where f = -j* , g being the distance from the rim of the trailing vortex sheet 
and h being the normal distance between two adjacent vortex sheets at the rim 
of the trailing vortex system. 

Prandtl's formula was derived for propellers with finite velocity of 
advance and negligible slipstream contraction. The basic concept, however, 
has been utilized to study the hovering rotor problem (Ref. 14) . Applying 
Eq. (D-l) to the ultimate wake, one has 

8 = R 1 " r l 
2TTU (R ) 

h = 

where N^ is the number of blades. Consequently, 

\ (1 - x) 

f = (D-2) 

The total circulation of the rotor is therefore 



and 
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^ TTCU l r l ^ 


(D-3) 


The thrust and the power are, in terms of r , 


T = pl^ 


-R 

% 

o 


F b r(n - §w) dr 


(D-4) 


p = pi^n 


>*R 

i 

o 


F, rudr 
b 


(D-5) 


Thus the thrust and power coefficients are 


c J 


- 8) 


ydy 


(d-6) 


■J 1 


rl^ r 


(D-7) 


where 


C = 


PrrftR* 


'(D-8) 


With the contraction factor k = , one has 


y = x 


u 

Hr 


2/~2~ 


and 


U) 1 , 

n " 2 b 


Thus, Eqs. (D-7) and (D-8) give 


38 



and 


T 

c 




yJy 


(d-9) 


p 

c 


— — J Caydy 
J~2 o 


(D-IO) 
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Figure 1 a General flow features 
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Figure 2(a) Optimum radial distribution of circulation 
for various values of thrust coefficient. 
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Figure 4 .- Comparison of approximate solution with present 

result - optimum angular velocity distribution 

for T = 0.010. 
c 

— ■ — — — approximate solution (Ref. 5 ) 

— — — — - ■ — present result 






Figure 6.- Comparison of approximate solution with present result 
optimum circulation distribution for T c = 0 .010 . 

— ~ “ ■— “ approximate solution (Ref. 5) 

. i ■ . . present result 




Figure 7.- Radial pressure distribution in the 
ultimate wake for T =0 .010 . 









